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Abstract 

We show that by means of connected-graph expansions one can effectively 
generate exact high-order series expansions which are informative of low-lying 
excited states for quantum many-body systems defined on a lattice. In par- 
ticular, the Fourier series coefficients of elementary excitation spectra are 
directly obtained. The numerical calculations involved are straightforward 
extensions of those which have already been used to calculate series expan- 
sions for ground-state correlations and T = susceptibilities in a wide variety 
of models. As a test, we have reproduced the known elementary excitation 

spectrum of the transverse-field Ising chain in its disordered phase. 
PACS numbers: 71.10.+x,05.50.+q,67.40.Db 
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Several years ago, Singh, Huse and this author described an efficient, general method 
for calculating high-order perturbations expansions for ground-state correlations and T = 
susceptibilities in quantum many-body systems defined on a lattice. (It is worth noting, 
belatedly, that we rediscovered methods that had been developed first by Hamer, Nickel, 
and others; a convenient entry point into that literature is He et al. This method had 
two principle facets: (1) A connected-cluster (a.k.a. "linked-cluster" , "connected-graph", or 
simply "cluster") expansion for the properties of interest, combined with (2) a formulation of 
Rayleigh-Schrodinger perturbation theory in terms of a recurrence relation, for the purpose 
of constructing the ground state energy and eigenvector perturbatively on each cluster. 
This method has been applied to a wide variety of problems involving S = 1/2 and 5* = 1 
antiferromagnets j| and also to the Kondo lattice model of mobile fermions and stationary 
spins ||. In every case one considers a family of Hamiltonians parameterized by A, 

H = H + XH U (1) 

in which the "unperturbed" Hamiltonian is trivially diagonalizable and of fixed degeneracy 
on any connected cluster of sites. 

The perturbation expansion approach is comparable in many respects to finite-size di- 
agonalization studies. A major advantage of the former is that the series expansions do 
not suffer from finite-size and shape effects — they are exact for the infinite lattice. A ma- 
jor disadvantage has been that many important properties associated with excited states, 
such as elementary excitation spectra and their associated residues, were inaccessible. Here 
I show that these properties can be calculated in terms of perturbation expansions about 
trivial unperturbed Hamiltonians, by direct (but hardly obvious, a priori) extensions of the 
two crucial facets of the method for ground state properties which are listed above. The 
algorithms are simple in principle, and if computer code for generating the ground-state 
series expansions for some model already exists then it is straightforward to modify it so as 
to yield excited-state data. 

The key to the approach is to construct the effective Hamiltonian within a degenerate 
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(for H ) manifold of excited states for each finite cluster. The formulae which yield per- 
turbation expansions for the effective Hamiltonians are derived by considering a similarity 
transformation 

S^HS = H cS (2) 

where H is the Hamiltonian restricted to some finite cluster, and written explicitly in the 
basis of the eigenstates of H Q . For definiteness we assume there are L degenerate states 
out of M states total, and it is useful to order the eigenstates of H so that the first L are 
those of the manifold. We can choose S so as to not mix the states within the degenerate 
manifold, that is, we take it to be the identity within the upper left L x L block. The 
matrix S is to be constructed so that in H eS there are vanishing matrix elements between 
the degenerate manifold and the other states, that is, H cS breaks up onto an L x L block 
in the upper left, describing the interactions between the "dressed" states of the manifold, 
and an (M — L) x (M — L) block, describing all the other (irrelevant) physics of the cluster. 
Traditionally (for example, the Foldy-Wouthuysen transformation of the Dirac equation), S 
would be written as e and the Hermitian operator T written as a power series in A, the 
exponentials would be expanded, and the vanishing of the appropriate H eS matrix elements 
would lead to a set of equations for the terms in the series for T. Rather than expanding 
exponentials, we construct H eS simply by moving S* -1 to the other side and viewing S as a 
set of column vectors We now have an equation in which the vectors ip® with 1 < I < L 
are not coupled to those with L + 1 < / < M, and we can now entirely neglect the part of 
H eS that lies outside the L x L block. Then expand everything in sight in powers of A: 

^> = £AtyJ?\ (3) 
k 

H cS (l',l) = Y,^ k Hf(l'J) (4) 

k 

and collect the terms proportional to A fc in column / to yield 

Ho4 ] + = EE 4 ] Hf-Ai', 0- (5) 

k'=0l'=l 
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From this equation one can now almost immediately write coupled recurrence relations 
for the terms in the power series expansions of H cS and the tp( l >. It is only necessary to 
impose a "normalization convention" (which affects the form of the final equations but not 
any of the physical results): it is most convenient to require for states \l) in the degenerate 
manifold of H that 

(ItyP) = 4,oV (6) 

Then one has for the effective Hamiltonian 

Hf(m,l) = (m\H 1 \^_ 1 ), (7) 

and for the projection of \if)"' > onto an eigenstate state of H outside the manifold the 
result is 

Hrf } > = tt^t (WilA) - E E Hf_ kl (l'J)(n\^; ] )) , (8) 
- V fc , =1 v=l J 

where Eq is the unperturbed energy of the states in the degenerate manifold and E n = 

(n\H \n). These last two equations are the desired recurrence relations. 

The next problem to address is how to carry out subgraph subtraction, so that one 
can obtain an exact nth order expansion by considering only a finite number of connected 
clusters. (See for example Ref. for a general discussion of subgraph subtraction, which is 
at the heart of the cluster expansion approach to constructing series expansions.) But at the 
moment, it is not entirely obvious what quantities one is supposed to be subtracting! The 
point one has to establish is that weight of the quantity Q under consideration vanishes for 
a disconnected graph; or, equivalently, for a graph C consisting of two unconnected clusters 
A and B, Q c = Qa + Qb- 

Since we are going to have to discuss specific quantities it is advantageous to consider a 
specific model, and for the sake of clarity we choose a particularly simple one: the S =1/2 
transverse-field Ising model defined by 

# = -E<-^E<^- ( 9 ) 

* (ij) 
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We will consider a "disordered state" expansion, in which H is the first term and Hi the 
second term (except for the factor A). Although the explicit discussion in the next several 
paragraphs will be concerned entirely with this model, it should be fairly clear that other 
models (such as the Kondo lattice model, variants of the t-J model, the Heisenberg-Ising 
model, and other models for which the ground-state properties have been studied by means 
of series expansions) may be treated along similar lines. 

For any iV-site cluster the ground state of H is the unique state with a* ot = N, namely, 
| | ... (AT times) |). The lowest excited states form an TV-fold degenerate set, the most 
natural basis set being | J,f . . . (N — 1 times) |), etc. One can label these "single-flipped- 
spin" states by the site of the cluster which carries the flipped spin. In the thermodynamic 
limit, one could equally well say that H has a dispersionless mode with energy 2. We will 
now proceed to show how one can use perturbation theory to determine how this spectrum 
evolves as the Ising coupling is turned on. 

Suppose we have a cluster C with two disconnected components A and B, as described 
above. (We will now attach A, B, or C as subscripts to a variety of quantities, to indicate 
that they are associated with one of these clusters). Consider the structure of H^ for the 
single-flipped-spin states. At order A , the states on which acts have a flipped spin in 
one of either A or B and no flipped spin in the other component, and it should be clear that 
even as A is increased from zero, and the states develop more structure, that there is no way 
for the spin-flip to "jump" between disconnected parts of the cluster. In other words, we 
assert that if we start with an unperturbed state in which there is a spin-flip in subcluster A 
but no spin-flip in B, and apply Hi (more precisely, the restriction of Hi to C) any number 
of times, the resulting state will have zero overlap with any unperturbed state in which there 
is a spin-flip in B but none in A. What this implies is that H^ has a block-diagonal form, 
specifically, 

Hf = [H cS + e B I] A © [H cS + e A I] B (10) 

where e A denotes the ground state energy of cluster A and / is the identity operator. Thus 
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H itself does not have a cluster expansion, since it is not simply additive for a disconnected 
cluster. However, if we subtract ec = e A + &b from all the diagonal elements on both sides 
of the equation above, 

[H eS - e c I} c = [H cS - e A I} A © [H cS - e B I] B , (11) 

we see that H eS — el does have a cluster expansion. In carrying out subgraph subtraction, 
one must keep track of which sites in the graph correspond to which sites in the subgraph, 
and the subtract all of the appropriate matrix elements, order by order. As a check of the 
validity of the calculation, one should find that the weight for H eS — el on a cluster with p 
terms of Hi should vanish identically up to order A p_1 . 

Once the weights have been evaluated for all graphs contributing up to a desired order 
in A, the excited state dispersion is trivial to evaluate. Each matrix element in the weight 
of H cS — el for each cluster is associated with some vector on the lattice. One must sim- 
ply sum all of the matrix elements associated with any given lattice vector r s ; denote the 
resulting sum of (power series of) matrix elements by t s . The effective Hamiltonian in the 
single-flipped- spin sector for the infinite lattice is now known to the desired order in A, and 
it is diagonalized by plane-wave eigenfunctions. The resulting dispersion relation for the 
excitation energy associated with a spin-flip is J2 S ts exp(iq • r a ). 

Just as in the case of ground-state expansions, one should take advantage of geometric 
and topological classification of the clusters in order to avoid redundant calculations: but 
that is a well-understood technical detail that requires no further discussion here. 

We have in fact implemented the algorithm described above for the S = 1/2 transverse- 
field Ising chain, and reproduced the exactly-known dispersion [1 + 2 A cos q + A 2 ] 1 / 2 to 
order A 11 (and for this model it would be easy to push the calculations further). Upon being 
informed of the algorithm, Singh J7J promptly carried out the analogous calculation for the 
single-flipped- spin excitations of the S = 1/2 square- lattice Heisenberg-Ising antiferromag- 
net; by appropriate extrapolation to the Heisenberg limit he directly obtains a value of the 
spin-wave velocity which is consistent with the current best estimates. 
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Similar calculations of elementary excitation spectra for other models would be straight- 
forward; however, they do not exhaust the potential applications of this method. The 
"eigenvectors" \ip^ 1 '} contain information on correlations in the vicinity of the excitation. 
The set of overlaps (in the context of the transverse-field Ising model, for concreteness) 
{il)\af,\ipW) , where \ip) is the perturbatively-constructed ground state of the cluster, should 
contain all the information needed to evaluate the quasiparticle residues associated with the 
elementary excitations. In fact a calculation of the residues has already been accomplished 
for the Heisenberg-Ising antiferromagnet, and will be reported separately. One subtle point 
worth emphasizing here is that the method of constructing H eS is a similarity transformation 
but not a unitary transformation. Insofar as the calculation of the spectrum is concerned, 
the distinction is irrelevant. However, the states 1^ ) are not orthogonal, and this fact must 
be carefully taken into account whenever they are used in further calculations. 

Another interesting extension would be to "two-particle" excited states. Here the sit- 
uation is rather more complicated than for "single-particle" excitations such as the one- 
flipped-spin states of the transverse field Ising model, on two counts. First, the two-particle 
effective Hamiltonian for the infinite lattice is not, in general, trivially diagonalizable like 
in the one-particle case. So the value of this method would be to yield a finite number of 
terms in an exact series expansion for the quasiparticle interactions: the next step would 
be along the lines of traditional many-body theory, or perhaps even exact diagonalization! 
(One would have dramatically reduced the size of the Hilbert space, compared to exactly 
diagonalization of the original model, but at the cost of introducing more complicated, long- 
range interactions.) Second, the analog of Eq. (|TTD is considerably more difficult to write 
down, and it is not clear precisely what quantity exhibits a cluster expansion. We believe 
these issues merit examination in future work, as well. 

At this point it is possible to sensibly compare the method outlined here with the few 
other series expansions methods for excited states which exist in the literature. Expansions 
for the elementary excitation spectra of the transverse-field Ising model on various lattices, 
and of the dimerized Heisenberg antiferromagnetic chain, were carried out some years ago 
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by Pfeuty and Elliott || and Harris ||, respectively. Neither of these methods is in the 
form of an explicit cluster expansion, nor do they appear to be readily implementable as 
a computer program, as required for high-order calculations. So far as we are aware these 
methods have not seen further use. 

Barber and Duxbury [[IT]] constructed the expansion for the lowest excitation energy in 
the next-neighbor transverse field Ising chain via Rayleigh-Schrodinger perturbation theory 
on finite rings. Because this is not a cluster expansion method it cannot be efficiently 
extended to higher-dimensional lattices, and it is not clear whether it can be adapted so as 
to yield the entire elementary excitation spectrum, either. However, this approach has been 
applied to other one- dimensional models (see for example [|TT|j ). 



Finally, there is a method which is due to B. G. Nickel [|12| and which has been applied 
by Hamer and coworkers to several spin models 0|13j- This method is almost a cluster 
expansion; to be precise a finite set of disconnected clusters need to be included as well as 
the connected clusters. It is also limited to the lowest excitation energy only. 

What sets the method presented in this paper apart from those already reported in the 
literature is the recognition that there is an excited state quantity which is additive over 
disconnected clusters, so that the usual machinery of connected-cluster expansions can be 
applied. Because this quantity turns out to be the effective Hamiltonian in the single-particle 
excitation subspace, one can obtain the entire elementary excitation spectrum as readily as 
the gap. 
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